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ABSTRACT. In this paper, a comparison of the Adomian and Taylor polynomial solutions of 
some nonlinear Ordinary Differential Equations of constant coefficients to those equations 
with variable coefficients is presented. The total derivatives of the nonlinear functions 
involved in the problem considered were derived in order to obtain the Adomian polynomials 
for the problems. For Taylor polynomials, the nonlinear functions involved were iteratively 
differentiated. Numerical experiments show that Adomian Decomposition Method can be 
extended as alternative way for finding numerical solutions to ordinary differential equations 
of variable coefficients. Furthermore, the methods are easy with no assumption and they 
produce accurate results when compared with other methods in literature. Moreover, the 
Taylor polynomial solution gives more accurate results compared to the results by Adomian 


polynomial solution in terms of the absolute errors. 


1.0 Introduction 
Differential equations generally and nonlinear differential equations in particular often do 
not have exact or closed form solutions, Batiha et al (2008). Nonlinear differential equations 
abound in many branches of applied mathematics such as psychology, chaos, potential theory, 
growth rate in biology and elasticity, Jimoh and Taiwo (2015). The Adomian Decomposition 
Method (ADM) proposed by G. Adomian in the 1980’s has been a subject of discussion 
and investigation by researchers, Almazmumy et al (2012) and Luo (2005). The (ADM) has 
been used by several mathematicians to solve nonlinear problems due to its simplicity and its 
ability to circumvent the problems of linearization and perturbation, Zhang et al (2006). The 
Adomian Decomposition Method involves separating the problem under consideration into 
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linear and nonlinear components, Hendi et al (2012). The nonlinear portion of the problem 
is decomposed into a series of polynomials. These are then called Adomian polynomials. 

Adomian decomposition method is a semi-analytical method for solving ordinary and par- 
tial nonlinear differential equations, Wazwaz (2000). The method, developed from 1970s to 
1990s by George Adomian is also applied to solve both linear and nonlinear Boundary Value 
Problems (BVPs) and integral equations. The numerical result is obtained with minimum 
amount of computation or mathematics, Hosseini (2006). Adomian technique is based on a 
decomposition of a solution of nonlinear functional equation in a series of functions. Each 
term of the series is obtained from polynomial generated by a power series expansion of an 
analytic function, Hosseini and Nasabzadeh (2007). Some of the advantages of the Adomian 
decomposition are that; it can be applied directly for all types of functional equations both 
linear and nonlinear and it has ability of greatly reducing the size of computational work 
while still maintaining high accuracy of the numerical solution, Hosseini and Nasabzadeh 
(2006). Adomian Decomposition Method (ADM) provides an analytical approximate so- 
lution for nonlinear functional equations in terms of a rapidly converging series, without 
linearization, perturbation or discretization, Hasan and Zhu (2009). 


Taylor series method is one of the singlestep methods which have different increment func- 
tions. From the application point of view, the Taylor series method has a major disadvantage. 
The method requires evaluation of partial derivatives of higher orders manually, Hoffman 
(2001). However, with the advent of computers, this challenge is being surmounted. The 
Taylor series can be written as the finite Taylor series, also known as the Taylor formula 
or Taylor polynomial. The Taylor series is a polynomial of infinite order, Vygodsky (1975). 
Thus, 


1 
f(x) = f (z0) + f'(z0)(£ — £0) + HF" (20) — 7 ee 
It is, of course, impractical to evaluate an infinite number of terms, Hoffman (2001). The 
Taylor polynomial of degree n is defined by 


f(x) = Pale) + Rnyi(*) 
where the Taylor polynomial P,,(x), and the remainder term R,,4;(x) are given by 
1 1 
P,(x) = f (£0) + f'(£0)(£ — £o) + ait (to) (2 = 29)" tant Sf (20) (a — zo)” 


1 
(n+ 1)! 
The Taylor polynomial is a truncated Taylor series, with an explicit remainder, or error term, 


Rn (£) = FOV (O\(@ — ao), Ea 





Mileties and Molnarka (2004). The Taylor polynomial cannot be used as an approximating 
function for discrete data because the derivatives required in the coefficients cannot be de- 
termined. It does have great significance, however, for polynomial approximation, because 
it has an explicit error term, Mileties and Molnarka (2005). 
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2.0 Adomian Polynomials 
Consider a functional equation 


(1) u=f+L(u) + N(u) 


where L and N are respectively linear and nonlinear operators and f is a known function. 
By Adomian decomposition method, the solution u(x,t) of (1) is decomposed in the form of 


an infinite series 
(2) Ue) = S| tn(2, t) 
n=0 


Furthermore, the nonlinear function N(u) assumes the following representation: 
(3) N(u) = Y An(tio, th, Un) 
n=0 


where, A’ s are nth order Adomian polynomials. In the linear case, N (u) = u then An simply 
reduces to tn. 

Cherruault and Adomian (1993) gave a method for determining these polynomials by parametriz- 
ing u(x,t) as 


(4) ual t,t) = So un(2, ta” 
n=0 
and assuming N (ua) to be analytic in a, which decomposes as 


(5) N (ua) = X. An (uo, tr, «++ Un Ja” 


n=0 
Hence, the Adomian polynomials A, are given by 


_ 10"N(ua) 


(6) An (uo, U1, JH in) = al Jar 


la-o0Vn € No 


where, Nm = [n € Nu0:n > m] and N denotes the set of positive integers. 
Rach (1984), suggested the following formulae for determining Adomian polynomials, Rach 
(2008): 


(7) Ao(uo) = N (uo) 


(8) An (Uo, U1, -3 Un) = c(k, n) NE (un )Yn € N 


k=1 
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The first-few Adomian polynomials are listed by Zhu et al (2005) as follows: 


Ao = f(t, Yo) 
A, = yı f'(t, yo) 
Az = yof'(t, yo) + surf" (t, Yo) 
As = ys f'(t, yo) + y1y2f" (t, yo) + zy? f" (t, Yo) 


where primes denote the partial derivatives with respect to y. 


(9) 


3 Adomian Polynomial Solutions of Ordinary Differential Equations 
The generalized first order nonlinear equation considered is given by 


(10) y' = f(x,y) 
together with the initial condition 
(11) Y(Ta) = Ya 


The main goal of this article is to extend the Adomian decomposition method, with modi- 
fication in order to obtain a polynomial solution of (10) and (11). Adomian Decomposition 
Method (ADM) solves nonlinear operator equations for any analytic nonlinearity, providing 
an easily computable, readily verifiable and rapidly convergent sequence of analytic approx- 
imate solutions. Since it was first presented in the 1980’s, Adomian decomposition method 
has led to several modifications on the method made by various researchers in an attempt to 
improve the accuracy or expand the application of the original method, Duan (2011). The 
choice of decomposition is non-unique and provides a valuable advantage to the analyst, per- 
mitting the freedom to design modified recursion schemes for ease of computation in realistic 
systems, Duan (2011). In order to obtain the Adomian polynomial solution of (10) and (11), 
the nonlinear variable coefficient equation (10) is written in its operator form as 


(12) Ly + Ry+Ny=F 


where F is a known function and y is the unknown function to be determined, L is the 
linear operator to be inverted, R is the linear remainder operator and N is the nonlinear 
operator, which is assumed to be analytic. The choice for L and its pair L7! (inverse of L) 
are determined by the equation being considered, hence the choice is non-unique. Here, we 


L is chosen to be 
d 


= 50) 
and thus its inverse L~! follows as the one-fold definite integration operator from zo to g. 
Thus, 


LLy=y- yY 
assumes the initial value as 


P = Ya 
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Remark: For nth-order differential equation, the choice of L is given by 
d” 
==) 
x 
and its inverse, L~' is the n-fold definite integration operator from xp to x. Thus, y absorbs 
the initial value as 


n-1 (x ee jk 
— To 
a ie 
k=0 
Applying the inverse linear operator L~! to both sides of equation (12), leads to 
(13) y = B(x) — L [Ry + Ny] 


where B(x) = y + LF 
The unknown function, y is expressed in a series of the form 


(14) y= Xo 


and the nonlinear term Ny is decomposed into a series 
k=0 


where the As which depend on yo,¥1,...,yx, are called the Adomian polynomials, and are 
obtained for the nonlinearity Ny = f(y) by 


k 
Sey hab ote 


(16) Ak = Toye 


where A is a formal parameter. 
However, in this work, the formulae (9) are modified to obtain expressions for the first-few 
Adomian polynomials Ap through Az, inclusively, as 


Ao = f (to, yo) 
A; = yı f’ (to, Yo) 
(17) Az = yof' (to, yo) + sy f" (to, Yo) 
Az = ysf'(to, yo) + yryot” (to, yo) + ey? f" (to, Yo) 
A4 = ya f'(to, yo) + $42? f" (to, yo) + Amys f” (to, yo) + Aui“ fO (to, yo) 
where primes denote total derivatives of f(t, y) at (to, yo). 
Using the Ais in (13) - (15), the recursive formula for yp41 is obtained as 


(18) T = AO E E E 
0 


4 Taylor Polynomial Solutions of Ordinary Differential Equations 

Barrio (2005), used a generalized Taylor series method for solving nonlinear ordinary differ- 
ential equations of constant coefficients to find a general expansion for a given function f(t). 
For Taylor polynomial solutions of first order ordinary differential equations, the initial value 
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problem given by an implicit differential equation (10) and initial value given by equation 
(11) is considered. 

Thus, if the solution g(x) of (10) and (11) is n times differentiable, the n-th degree Taylor 
polynomial of g(x) is then sought for. A typical situation is an initial value problem (10) 
and (11) for a complex function g with analytic right hand side f. In this case, the solution 
g is analytic in a neighbourhood of the initial point 79. The Taylor coefficients are given by 


n lk) (x 
(19) Tn(g, £, £0) =+? I 0) 
k=1 ` 





(£ — zo)" 


with the derivative list as 


g(x) = f(@,Wly=ote) = fey) 


ð ð 
g(x) = so, y) + L Caya) 
ð ð 
Z Îl a, y) + Lia y) f(z, Y) ly=ote) = fo(a, y) 





(20) gf) = F(a, y) + PRA la, y) fle, vlys = Falay) 
In the case of an implicit second order differential equation 

(21) y" = f(z,y, y’) 

with initial conditions 

(22) Y(Ta) = Ya 

and 

(23) y (a) = W 


the n-th degree Taylor polynomial solution has the form 


g™® (xo) 
k! 





(24) Tn(g,£, £0) = go + gı(x — xo) + ` (x — x9)* 


k=2 


and the same method as before after calculation of g) (29) is applied for k > 2. 
If the right hand side of (21) is a function of the three variables x, y and u, then the 
derivatives of g(x) are obtained iteratively as 


g(x) = f(x, g(x), g'(x)) 
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and by the chain rule 


o o o 
(2) = Elen + Ale yuge) + Ea, yuga) 


20h of Of 
7 gr ho YU) T dy (x,y, uju + gu $ V WFE, y, Wlu=g'(a) 


for y = g(x) and iteratively 


o fk—ı o fk-1 O fr- 
(k) — ý ana E 
g (x) Ox (x,y, U) a dy (x,y, u)g (x) a du Foe) ay (x) Fy(a,y, u) 











and g% (x) is evaluated at the point x = 2 to get g% (£o) by taking the limit u — yı, 
Y — yo and x —+ zo which yields the result. 

After obtaining the polynomials g“) (xo) for k = 1, 2, 3,...,n, they are then substituted into 
equation (19) for the Taylor polynomial solution of equations (10) and (11). 


5 Evaluation of the Error 
In this paper, error is defined as 


Error = max |ExactValue — ApproximateV alue| 


In case the exact solution is not available, the results are compared with those in literature. 


6 Illustrative Examples 

The Adomian and Taylor polynomial solutions obtained for some nonlinear first order ordi- 
nary differential equations of variable coefficients are presented in this section. The results 
obtained are tabulated for comparison. 


Example 1: Consider the first order nonlinear differential equation 


(25) y Say +1 
with the initial condition 
y(0) = 1 
(Griffths and Higham, 2010) 
Example 2: Consider the first order nonlinear differential equation 


(26) y = ry’ 


with the initial condition 


The exact solution is 
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Solution by the Adomian Decomposition method: 
y(x) = 0.98973 + 0.1693672 — 1.64583x" + 13.43752° — 41.66672* + 52.08332° 


(Griffths and Higham, 2010) 
Example 3: Consider the first order nonlinear ordinary differential equation 


(27) y =r +y 
with the initial condition 
y(0)=1 


Solution by the Adomian Decomposition method: 





(Jain et, al, 2012) 


Example 4: Consider the first order differential equation 


yY 5 2 3 
28 Pen Bi 
(28) y =z 73T 
with the initial condition i 
1 = —— 
y(1) Ti 


Solution by the Adomian Decomposition method: 
y(x) = —1.34642 + 7.334842 — 11.4929x + 10.6446x° — 5.44012* + 1.300002° 


(Griffths and Higham, 2010) Tables of Results 
TABLE 1 
Numerical Results for Example 1 






































x | Literature | Adomian Poly. | Taylor’s Poly. 
0.0) 1.00000 | 1.00000 1.00000 
0.1 | 1.10572 | 1.10593 1.10569 
0.2 | 1.22600 | 1.22832 1.22579 
0.3 | 1.386856 | 1.37617 1.36543 
0.4} 1.54210 | 1.56171 1.53077 
0.5) 1.76596 | 1.80078 1.72917 
0.6 | 2.06740 | 2.11336 1.96936 
0.7 | 2.50470 | 2.52396 2.26170 
0.8 | 3.20741 | 3.06208 2.61834 
0.9) 4.53927 | 3.76270 3.05344 
1.0 | 8.00782 | 4.66667 3.58333 




















Literature: ‘Taylor Series Method 
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TABLE 2 
Numerical Results for Example 2 































































































x Exact Solution | Adomian Poly. Error Taylor’s Poly. Error 
0.10 1.00000 1.00000 0.0000E-5 | 1.00000 0.0000E-5 
0.11 1.00105 1.00106 1.0000E-5 | 1.00105 0.0000E-5 
0.12 1.00220 1.00223 3.0000E-5 | 1.00221 1.0000E-5 
0.13 1.00346 1.00348 2.0000E-5 | 1.00346 0.0000E-5 
0.14 1.00482 1.00485 3.0000E-5 | 1.00482 0.0000E-5 
0.15 1.00629 1.00631 2.0000E-5 | 1.00629 0.0000E-5 
0.16 1.00786 1.00789 3.0000E-5 | 1.00786 0.0000E-5 
0.17 1.00954 1.00956 2.0000E-5 | 1.00954 0.0000E-5 
0.18 1.01133 1.01136 3.0000E-5 | 1.01133 0.0000E-5 
0.19 1.01322 1.01324 2.0000E-5 | 1.01322 0.0000E-5 
0.20 1.01523 1.01527 4.0000E-5 | 1.01523 0.0000E-5 
TABLE 3 

Numerical Results for Example 3 

x Lit. Adomian Poly. | Taylor’s Poly. 

0.0 | 1.000000 | 1.00000 1.00000 

0.1 - 1.11151 1.11146 

0.2) 1.25253 | 1.25367 1.25292 

0.3 - 1.44209 1.43837 

0.4 | 1.69318 | 1.69892 1.68749 

0.5 - 2.05417 2.02708 

0.6 | 2.62293 | 2.54694 2.49251 

0.7 - 3.22677 3.12913 

0.8 | 5.49363 | 4.15486 3.99375 

0.9 - 5.40536 5.15604 

1.0 | 34.3075 | 7.06667 6.70000 




















Lit. Taylor series method of order four 





10 A. K. JIMOH, A. M. AYINDE 


TABLE 4 
Numerical Results for Example 4 





















































x Lit. | Adomian Poly. | Taylor’s Poly. 
1.0 | 1.00000 | 1.00000 1.00000 
1.1 - 1.11227 1.11145 
1.2 | 1.25478 | 1.25371 1.25289 
1.3 - 1.44140 1.43833 
1.4 | 1.69912 | 1.69808 1.68744 
1.5 - 2.05371 2.02702 
1.6 | 2.57034 | 2.54703 2.49243 
1.7 - 3.22713 3.12903 
1.8 | 4.27532 | 4.15499 3.99363 
1.9 - 5.40508 5.15590 
2.0 | 7.20134 | 7.06686 6.69983 
Lit. Modified Euler Method 





Graph of Example 1 
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Figure 1: The behaviour of the Taylor Series Method compared with 
the Adomian and Taylor polynomial solutions. 
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Graph of Example 2 
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Figure 2: The behaviour of the exact solution compared with 
the Adomian and Taylor polynomial solutions. 


Graph of Example 3 
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Figure 3: The behaviour of the Taylor Series Method compared with 
the Adomian and Taylor polynomial solutions. 
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Graph of Example 4 
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Figure 4: The behaviour of the Modified Euler’s Method compared with 
the Adomian and Taylor polynomial solutions. 


7 Discussion of Results 
The Adomian and Taylor polynomial solutions of some nonlinear ordinary differential equa- 
tions of variable coefficients were obtained. The results were compared with the exact so- 
lutions (where available) and some existing results in literature. The absolute errors show 
that the results by the two polynomial solutions are in excellent agreement with the exact 
solutions. Moreover, the Taylor polynomial solution gives a better approximation than the 
Adomian polynomial solution as can be seen in the tables 1 - 4 presented. It is also noted 
that the Taylor polynomial solution performs better than the Taylor series method as the 
results by the Taylor series method diverge as the value of x moves away from the initial 


point. 
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